Uncertainty principle with quantum Fisher information * 



Attila Andai f 
RIKEN, BSI, Amari Research Unit 
2-1, Hirosawa, Wako, Saitama 351-0198, Japan. 

October 10, 2007 



Abstract 

In this paper we prove a nontrivial lower bound for the determinant of the covariance matrix 
of quantum mechanical observables, which was conjectured by Gibilisco, Isola and Imparato. The 
lower bound is given in terms of the commutator of the state and the observables and their scalar 
product, which is generated by an arbitrary symmetric operator monotone function. 

Introduction 

The basic object in the statistical description of a classical physical system is a probability space 
(Q,,B,n), where the measure \i determines the state of the system and the physical quantities are 
measurable fi — ► R functions. The covariance of the quantities X, Y € L 2 (il,M., /i) is defined as 



Cov^(X,Y) = f XYdH - ( I XdA ( [ Y dfJ. 

Jo, \Jn / \Jn 



and the variance of a quantity is Var^(X) = Cov At (X, X). The Cauchy-Schwartz inequality in this 
setting gives 

Var^pOVar^Y) - Cov^F) 2 > 

which can be reformulated as 

fCoy^X,X) Cav^X,Y) 
ael \Cov^Y,X) Cov (I (7T) 

The quantum mechanical Hilbert space formalism gives a mathematical description of particles 
with spin of Concentrating on the spin part of non relativistic particles one can build a 

proper mathematical model in an n dimensional complex Hilbert space. This is the simplest physical 
realization of an n-level quantum system. The states of an n-level system are identified with the set of 
positive semidefinite self-adjoint n x n matrices of trace 1, and the physical observables are identified 
with the set of self-adjoint n x n matrices. For a given state D the (symmetrized) covariance of the 
observables A and B is defined as 

Cov D (A, B) = - (Tr(DAB) + Tr(DBA)) - Tr (DA) Ti(DB) 
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and the variance of a observable is Varu(A) = Covd(A, A). From the Cauchy-Schwartz inequality 
we have 

Vat D (A)Vax D (B) -Cov D (A,B) 2 > i| Tr(p [A, B])\ 2 

which is known as the Schrodinger uncertainty principle [19) . Without the covariance part, one gets 
the Heisenberg uncertainty relation [5] . The Schrodinger uncertainty principle can be reformulated as 



det 



fCov D (A,A) Cov D (A,B) 
{ Cov D {B, A) Cov D (B,B] 



> det 



1 fTr(D [A, A]) Tr(D [A, B]) 

2 \Tr(D [B,A}) Tr(D[A,A]) 



This form was generalized by Robertson for the set of observables (Ai)i t „ n as 



detl [CovoiA^Aj)] 



h,j=l, 



,N 



> det 



In this formula the lower bound is given by the commutators of the observables. 
In this paper we prove the inequality 



det (•oY l> {A i ,.A l i h j , v > det 



M (i[jD)A])i [ AA . ])D/ 



h,j=l, 



,N 



(1) 



where the scalar product (•, -)dj is induced by an operator monotone function /, according to Petz 
classification theorem [18]. The inequality |T]) was studied first just in the case N = 1 for special 
functions /. The cases f(x) = Jsld{x) = i±£ anc | f( x ^ = f WY (x) = \{\/x + l) 2 were proved by 
Luo in [Hjrrj]. The general case of the Conjecture was proved by Hansen in [7] and shortly after by 
Gibilisco, Imparato and Isola with a different technique in [3]. 

In the case N = 2 the inequality was proved for / = fwY by Luo, Q. Zhang and Z. 
Zhang [T3J EH EHj- The case of Wigner-Yanase-Dyson metric, where fpix) = r§B^W^W~i) 
(0 G [—1, 2] \ {0, 1}) was proved independently by Kosaki [10] and by Yanagi, Furuichi and Kuriyama 
|20j . The general case is due to Gibilisco, Imparato and Isola p2[B]. Gibilisco and Isola emphasized 
the geometric aspects of the inequality |T]) and conjectured it for general quantum Fisher information 

In a recent paper Gibilisco, Imparato and Isola proved the inequality in the N — 3 real case [1] 
and conjectured that 



det 



/(0) , 



i[D,A h ],i[D,Aj\) 



Jl/DJ 



h,j=l, 



,N 



> det 



hj=l, 



.N 



(2) 



holds too in general if jj^y > They proved this conjecture for every N £ N and for every 

appropirate function / recently in [5J. 

In this paper we prove this inequality too in a bit stronger form. 



1 Quantum Fisher information 

The states of an n-level system are identified with the set of positive semidefinite self-adjoint n x n 
matrices of trace 1. The states form a closed convex set in the space of matrices and denote by M n 
its interior, the set of all strictly positive self-adjoint matrices of trace 1. Let M^°i a be the real vector 
space of all self-adjoint n x n matrices of trace 0. 

The space M. n can be endowed with a differentiable structure [5]. The tangent space To at 
D € M n can be identified with M n °l a . A map 

K : M n x Mi% x Mi% - C (D, X, Y) i-> K D (X, Y) 
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will be called a Riemannian metric if the following condition hold. For all D G M. n the map 

K D : Mi°l x Mi% - C (X, Y) -> D (X, Y) 

is a scalar product and for all X G M^°i a the map 

K.{X,X) : M„ -> C Dh^if D (X,X) 

is smooth. 

Let (lf( m ')fneN be a family of metrics, such that is a Riemannian metric on A4 n for all m. 

This family of metrics defined to be monotone if 

K { ^ D) {T{X),T{X))<K { £\X,X) 

for every stochastic mapping, that is completely positive, trace preserving linear map T : M n (C) — > 
M m (C), for every D G .M„ and for all X G M^°i a and for all m, neN. 

Theorem 1. Petz classification theorem '181. There exists a bijective correspondence between the 
monotone family of metrics (K^) ne jq and operator monotone f : R + — > K functions such that 
f(x) = x/(x _1 ) hold for all positive x. The metric is given by 

K%\X,Y) = Tr(x(Rl D f(L ntD R-^)Rl D )-\Y)) (3) 

for all n G N where L„ iD {X) = DX , R n . D (X) = XD for all D,X G M„(C). 

For simplicity we denote by (•, - ) D ^ the scalar product given by Equation ([3]) in the tangent 
space of the point D. Let denote by F op the set of operator monotone functions / : R + — * R with 
the property f(x) = xf(x~ 1 ) for every positive parameter x and with the normalization condition 
/(l) = 1. Here are some elements of the set J- op from Refs. [151 H7] : 

1 + x 2x x-1 2(x-l) 2 2(x-l)s/x 2x a+1 ' 2 [3(1 - (3){x - 1) 2 
2 ' 1 + logx ' (1 + a;) (log a;) 2 ' + logx' 1 + x 2a ' (x/ 3 - l)^ 1 -^ - 1)' 

where < a < 1/2 and < < 1. We also introduce the sets 

= {/ G I /(°) ^ °> and ^op 5 = {/ G ^"op I /(0) = 0} . 
Theorem 2. /5/ For / G ^ and x G R+ set 

/» = i( (l+ i)- (l -i)^). 

For a function / G .Fop we define 

m f (x,y) = x/ (|) 

which will be used in sequel. 
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2 Volume uncertainty principle for Fisher information 

The observables of an n-level quantum system are identified with the n x n self-adjoint matrices, 
which we will be denoted by M„ !Sa . For observables A,Be M„. sa , state D £ A4 n and function 
/ G J-™ p we define the covariance of A and B with respect to D and the quantum covariance of A and 
B with respect to D and / as 

Cov D (A,B) = -(Tr(DAB) +Tr(DBA)j - Tr(DA) Tr(DB) 

Qov A/ (A, B) = ffl (i [D, A] , i [D, B]) D } . 

Theorem 3. Let {e h ) h=1 be a complete orthonormal base composed of eigenvectors of D € M n , 
and (Xh) h= i n be the corresponding eigenvalues. For matrices A, B S Af„. sa we associate matrices 
A and *B whose entries are given respectively by *Ahj — (Aoe h , ej) , *Bhj = (Boe hl ej). We have the 
following identities. 



1 " 

Cov D (A, B) = ReTr(ZMo-Bo) - ^ £ ( Xh + A i) M* A hj%h) 

hj=l 

^ n n 

Qov Df (A 1 B) = - ^ (A/j + \j)Ke(*A h j*Bj h ) — ^ m f ~(X h , Xj) Re(* A hj*B jh ). 



2 

h,j=i h,j=i 

In the equations of the previous Theorem the real part function Re can be omitted, since the 
matrices *A and *B are self- adjoint 



hi 



*A jh = (A Q e h ,ej) = {e^Aoe^ = (Aoe^e^) = *A hj 
and for real, symmetric coefficients a^j 

n n n 

^ a hj A hj *B jh = 2J a hj Re(*A hj *B jh ) +i 2J a hj lm(*A hj *B jh ) 

h,j=l h,j=l hj = l 

n n 

= ^ a hJ Re(*A h j*Bjh) + i ahh lm -(* A hh*B hh ) + i ^ a hj (im^/S^) + ka(*A jh *B, 

h.j=l h=l l<h<j<n 

n 

= ^ a hj Re(*A hj *B jh ) +i J2 a U lra(*A hj *B jh + *Ahj*B jh ) 

h,j=l l<h<j<n 
n 

= J2 a hj Re(*A hj *B jh ) 

h,j=l 

holds. So we will use the equations 

Cov D {A,B) = £ f^y^) *A hj *B jh (4) 
h,j=i ^ ' 

Qov DJ (A, B) = J2 (*h±hL - m f (\ h , Xj)j *A hj *B jh (5) 

for covariances. 
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We will call the set of nonzero matrices (Ai) ie i offdiagonally independent if none of them can be 
written as a linear combination of the other matrices and a diagonal matrix. We call the set (j4j)j e j 
offdiagonally dependent if there are nonzero (a,)iei coefficients such that the matrix ^2 ieI a%Ai is 
diagonal. 

Theorem 4. Consider a state D e A4 n , functions f,g e such that 

ffi > ff(0) 
f(t) g(t) 

holds and an N-tuples of nonzero matrices (j4^)/c=i,...,iv € M„ jSa . Define the N x N matrices Covp 
and Qov D j with entries 

[Cov D ] y =Cov2,(AW,^W)) 
[Qov Di/ ]..=Qov Di/ (^W,^)). 

(1) FFe /iai>e 

det(Cov D ) > det(Qov D y ) + dct(Cov D - Qov D y ) + i?(£>, /, N), (6) 
det(Qov DJ ) > det(Qov D J + det(Qov D >f - Qov D J + R(D, f, g, N), (7) 

where 

1 ^z} /m\ k_ N-k 

R(D, f,N)=J2[ k ) [det(Qov Di/ )] N [det(Cov D - Qov DJ )] » , 

/, fl , ^V) = E ( fc J [det(Qov Diff )] " [det(Qov A/ - Qov^)] " . 

(2) The following conditions are equivalent. 

a. det(Cov£>) = det(Qov D j) 

b. det(Qov D / ) = det(Qov D g ) 

c. The set of matrices (A^)k=i,... t N are linearly dependent. 

Proof. The matrix Qov^, j is obviously real and symmetric. First we prove that Qov^, j.- is positive 
definite. Let us define for indices 1 < h, j < n 

(f) Ah + Aj (o) \h + Aj . . 

a ^ = 2 ~ fc ' j) a « = 2 " ( ^ 

Elementary calculations show that = and for different indices > 0. Consider a vector 
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x € C N and define an n x n matrix as C — J2a=i x a *A^ a \ Then we have 

N 

(x,Qpv Dtf x) = wQov D!/ (i (a U (6) ) (9) 
0,6=1 

N n 

= E E 

a.b—l h,j—l 

n / iV \~ / AT \ 

/»,j=l \a=l / \6=1 / 

n 

= E 4?i^-i 2 >o. 

h,j=l 

We note, that if the set of matrices (*A^)k=i,...,N are offdiagonally dependent then there exists 
a nonzero vector x £ C N , such that C — Yl a =i x a*A^ is a diagonal matrix. For this vector 
(x,Qov£, j x) = holds, and it means that is an eigenvalue of the matrix Qovp ^, therefore its 
determinant is zero. If the matrices (*54W)fc=i,...,jV are offdiagonally independent then for every 
nonzero vector x we have (x, Qov^, ji) > 0, and in this case det(Qov £) j) > 0. So we have the 
following equivalence: 

det(Qov D f ) = if and only if CA.^)fe=i,...,iv are offdiagonally dependent. (10) 

We can repeat our arguments from Equation §§§ for the matrix Covd — Qov D * using a' ' — cS^ 
instead of . This lead us to the conclusion that Cov£> — Qovp t is real, symmetric, positive definite 
matrix. Since — a^fj > we have the following equivalence: 

det(Covc — Qov£, A = Q if and only if (A^)k=i....,N are linearly dependent. (11) 

If for functions f,g£ Top 

no) . 9(0) 

fit) g(t) 

holds for every positive parameter t then we have mj(x,y) < m,g(x,y) for every positive x and y. It 

means that a[ f ) - a ( $ > for every indices 1 < h.j < n. Using the previous arguments we conclude 
the positivity of the matrix Qov D j — Qav D g and the equivalence: 

det(Qov £) j — Qov D ) = if and only if (Aq )fc=i,... 1 jv are linearly dependent. (12) 

Using the Minkowski determinant inequality (see for example [T] p. 70.) for real symmetric 
positive matrices Qov D j and (Cov^ — Qov D f) we have 

[det(Qov A/ +(Cov D - Qpv DJ ))] * > [det(Qov D>/ )] * + [det(Cov D - Qov D>/ )] * (13) 
and for matrices Qov^, and (Qovp j — Qov D ) we have 

[det(Qov As +(Qov A/ - Qov Dig ))] * > [det(Qov D!9 )] * + [det(Qov 0>/ - Qov D g )] * . (14) 
These equations implies the first part of the Theorem. 
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To prove the second part of the Theorem assume that det(Cov£>) = det(Qov0 A Using 

the inequality (J6j) we get that det(Cov£> — Qov D *) — 0, which implies that (A^)k=i n are 

linearly dependent according to the equivalence (fTTj) . If (.Aq )fc=i,...,./v are linearly dependent then 
det(Cov£>) = and det(Qov D j) — 0, so in this case det(Cov£i) = det(Qov D j). Now we proved the 
equivalence of the a. and c. statements. The equivalence of b. and c. can be proved similarly. □ 

Theorem 5. Consider a state D S M. n , functions /, <? € such that 

/(0) ^ 9(0) 
/(*) g(t) 

holds and an N -tuples of nonzero matrices (AV*')t—i ,...,jv € M„. sa . Then for every t £ [0, 1] parameter 
the inequalities 

det{t Cov D +(1 - 2t) Qov D f ) > (1 - t) N det(Qov Z5 f ) + t N det(Cov D - Q,ov D / ) + R{D, f, N, t), 

(15) 

det(t Qov D f +(1 - 2t) Qov D g )> (1 - <) w det(Qov D g ) + < w det(Qov D / - Qov D g ) + R(D, f, g, N, t), 

(16) 

hold, where 

R(D, f, N,t)=J2 ( N k ) ((1 " *) ^/det(Qov A/ )) * (i ^det(Covr> - Qov A/ )) 

k=l ^ ' 

R(D,f,g,N,t)=J^ (*f\ ((1 - i) ^det(Qov A9 )) fe (t ^/det(Qov A/ - Qav Afl ))"~*. 

Proof. The Theorem is just an application of a generalized form of the Minkowski inequality, which is 
due to Firey [2]: for k x k real, symmetric, positive definite matrices K, L and for parameter t S [0, 1] 

dct((l - t)K + tL)i > (l-t)det(K)i + tdet{L)i. 

□ 

The t = \ case in the previous Theorem gives back Theorem ^ . 
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